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Abstract 

Conventional cloaking based on Euclidean transformation optics requires that the 
speed of light should tend to infinity on the inner surface of the cloak. Non-Euclidean 
cloaking still needed media with superluminal propagation. Here we show by giving an 
example that this is no longer necessary. 



1 Introduction 



10 and metamaterials 11 16 had lead to the 
18 and invisibility cloaks |4l|5l. Since then 



The development of transformation optics 
theoretical prediction of invisible dipoles 
the principles behind invisibility have been demonstrated in a number of experiments that 
feature cloaking in some reduced form. Electromagnetic cloaking was first realised for mi- 
crowaves of one frequency and polarisation 



20-26 



19 



Numerous subsequent electromagnetic ex- 
A carpet cloak does not make things invisible, 
Approximate cloaking has also been realised through 
Invisibility by plasmonic covering 



periments achieved "carpet cloaking'' 
but makes them appear to be flat, 
tapered waveguides 

well |29|. Furthermore, non-electromagnetic forms of cloaking have been demonstrated for 
acoustic waves 



17 was demonstrated as 



30,31 



However, the realisation of electromagnetic cloaking suffers from a practical and a fundamental 
problem 32 34 . The practical problem is that the material requirements for electromagnetic 



cloaking are very difficult to meet, whereas the fundamental problem is that perfect invisi- 
bility j5] requires that light should propagate in certain cloaking regions with a superluminal 
phase velocity that tends to infinity. In principle, this can be achieved by metamaterials, but 
only for discrete frequencies that correspond to the resonant frequencies of the cloaking mate- 
rial [6)[34] 



In contrast, acoustic cloaking 30,31 is much easier to achieve, since acoustic waves 
are not effected by relativistic causality and thus they are not restricted by their velocity in 
air. 



Broadband cloaking based on non-Euclidean geometries has been proposed in 32 , 33 1 to avoid 



the requirement for infinite light velocities for electromagnetic cloaking. In this proposal the 
speed of light is finite in the entire cloaking region and, therefore, this cloak has the potential 
of working for a broad range of frequencies. However a fundamental problem still remains. 
Since space has to be expanded to make room for the invisible region within the cloak, the 
implementation of this device will still demand superluminal propagation (i.e. propagation 
with a velocity that exceeds the speed of light in vacuum). The same is true for "carpet 
cloaking" 20 26 where the velocity of light in the cloaking device must exceed the speed of 
light in the host material, i.e. vacuum if such devices were to find practical applications. 
In this paper we give an example of a device that achieves complete electromagnetic cloaking — 
not just "carpet cloaking" — while all light velocities within the cloak are finite and less than 
the speed of lightfj] Through this example we demonstrate that invisibility cloaking is possible 
without superluminal propagation and anomalous material requirements. 



2 Problem 

Since transformation optics establishes a one-to-one correspondence between spatial geome- 
tries and dielectric media, invisibility can be visualised in terms of virtual geometries. There- 



lr rhe preprint 35 proposes a different method for cloaking without superluminal propagation. 
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Figure 1: A Euclidean cloaking device expands a single point in virtual space (blue dot in A) 
into an extended region in physical space (blue circle in B) through the curved transformation 
of the virtual coordinate grid. A light ray that is smoothly guided around the invisible region 
in physical space (B) appears to have passed through empty space (A), making the region 
within the blue circle (B) invisible. However, along the blue circle the phase velocity of light 
goes to infinity. 



fore, we explain the main ideas of this paper in terms of pictures. The complete calculations 
behind these pictures can be found in the Appendix. 

First, let us contrast our method to conventional cloaking based on coordinate transforma- 
tion. Figure 1 explains the conventional cloaking method described in [5]. A single point in 
virtual space - also called electromagnetic or optical space - is expanded into a finite region in 
physical space. The expansion is effected by the curved transformation of the straight virtual 
coordinate lines. Consequently, a light ray that traces out a curved trajectory in physical space 
(Fig. IB) appears to propagate along a straight line in virtual space (red line in Fig. 1A). This 
creates the illusion that light has propagated through empty space, making the central region 
in physical space (region bounded by the blue circle in Fig. IB) invisible. 
However, there is a fundamental problem with this device. Light crosses the central point 
in virtual space (blue dot in Fig. 1A) in an infinitely short time. This point turns into an 
extended region in physical space, which still has to be traversed in an infinitely short instant. 
Therefore, the phase velocity of light will tend to infinity along the outer boundary of the 
invisible region in physical space. Infinite light velocities can be achieved by metamaterials, 



but only for discrete resonant frequencies of the material 32 34 , fundamentally restricting 
the applicability of this cloak. 

In contrast, the broadband non- Euclidean invisibility device proposed in [32] achieves full 
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cloaking without the need for infinite light speeds. It makes use of the rich structure of non- 
Euclidean geometries to avoid infinite expansions and singularities. 

The background geometry of this cloak is established through the expansion of space, which 
is similar to optical conformal mapping |4([9)[34)[36] (see Appendix [Of . A finite region m 
physical space (red eye-shaped region in Fig. 2A) is turned "inside out" to produce a hidden 
infinite virtual plane (red lower plane in Fig. 2B). The blue outer region in physical space 
(Fig. 2A) is stretched to make the upper sheet complete without the red region (Fig. 2B). 
The two sheets remain connected by the yellow branchcut (Fig. 2B), through which light can 
pass from one sheet to another. The two virtual sheets provide the background geometry of 
the device. By the appropriate modification of the lower sheet, objects placed there can be 
hidden from sight. For the time being we will focus our attention on the blue region of the 
upper sheet, which we shall call the "stretched region" of the cloak. 

It is clear from Fig. 2 that the blue region was stretched everywhere by merely finite factors. 
Therefore, establishing this geometry will not require infinite light speeds in contrast to the 
previous cloaking example (Fig. 1). However, a problem still remains. A light ray that prop- 
agates along the straight yellow branchcut in virtual space (Fig. 2B) will have to trace out 
a slightly longer curved path along the physical branchcut (Fig. 2A) in the same amount of 



time. Since light travels in virtual space with the speed of light in vacuum 32 , the phase 



velocity in physical space will have to be faster than c. This can be achieved only for certain 



frequency bands that are around the resonant frequencies of the material 46 , which restricts 
the applicability of the cloak. 



3 Solution 

3.1 Geometry 

In the previous section we saw that the Non-Euclidean Cloak requires that light should propa- 
gate faster than the speed of light in vacuum. In order to avoid such anomalous light velocities, 
we construct a device that makes use of the same space expansion as the Non-Euclidean Cloak 
(Fig. 2), but is built up from different constituents. 

To slow the light rays down in the stretched region of the cloak (blue region in Fig. 3A), we 
place an optically dense medium (that is, a medium with high refractive index) on the upper 
sheet of virtual space (Fig. 3B). This medium must itself be invisible to avoid revealing the 
presence of the invisibility device. Therefore, we use the partially transmuted Invisible Sphere 
as the background profile for the upper sheet of virtual space (purple circle in Fig. 3B). All 
rays entering this Sphere perform a complete loop in it and emerge with their original direc- 
tion restored (Fig. 3B). Therefore, the entire sphere appears invisible. At the same time, if 
a light ray crosses the yellow branchcut while propagating through the Sphere, it will pass 
onto the lower virtual sheet of the geometry (Fig. 3B). We place Maxwell's Fish Eye Lens on 
the lower virtual sheet and surround it by a circular mirror (dark cyan circle in Fig. 3B) to 
guide all light rays that enter the lower sheet back to the yellow branchcut. The light ray 
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Figure 2: The expansion of space. The red eye-shaped region in physical space (A) is expanded 
into a red infinite virtual plane (B). The blue region (A) is stretched to make the upper sheet 
complete without the red region (B). A light ray propagating with the speed of light along 
the straight yellow virtual branchcut (B) will have to follow the longer curved yellow line in 
physical space (A). Therefore, in physical space light will have to propagate faster than the 
speed of light in vacuum. 



reemerges from the branchcut on the upper sheet and finishes its trajectory in the Invisible 
Sphere (Fig. 3B). Since the original direction of the ray is restored when it emerges from 
the device, it appears to have propagated through empty space. Note that there is a region 
outside the mirror on the lower sheet that is inaccessible to light. Any object placed in that 
region is undetectable to light and thus becomes invisible. Fig. 3A shows the light trajectory 
in physical space. The light ray (red) completes a loop that is characteristic of the Invisible 
Sphere, bounces off the spindle-shaped mirror (dark cyan) twice and leaves the device with 
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Figure 3: Light propagation in physical space (A) and virtual space (B). We place the partially 
transmuted Invisible Sphere on the upper sheet of our virtual geometry (region bounded by the 
purple circle in B) and put Maxwell's Fish Eye Lens on the lower virtual sheet and surround 
it by a circular mirror (dark cyan circle in B). Light explores only the interior of the mirror 
in virtual space (B). The unexplored infinite region of the lower virtual sheet (B) transforms 
into a finite region in physical space (white spindle-shaped region in A) that is unaccessible to 
light and where objects can be hidden. In physical space light completes a loop in the device, 
bounces off the mirror twice and leaves the device with its original direction restored. 



its original direction restored. Note that the infinite region outside the circular mirror on the 
lower virtual plane (Fig. 3B) transforms into a finite spindle-shaped region in physical space 
(white region in Fig. 3A). This region is inaccessible to light and any object placed in there 
remains hidden from sight. Since the stretched region of the cloak now sits in the optically 
dense Invisible Sphere, all light rays within it are propagating slower than the speed of light. 
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3.2 Building blocks 



In this section we explain the partially transmuted Invisible Sphere and Maxwell's Fish Eye 
Lens with the mirror, which serve as the two fundamental building blocks of our virtual ge- 
ometry. 

The Invisible Sphere is an optical lens that is invisible in the limit of geometrical optics [9 42 



Each incident ray completes a loop in the Sphere and emerges in the direction of its entry 
(Fig. 4A), creating the illusion that space is empty. 




Figure 4: The Invisible Sphere and its partial transmutation. Each ray entering the Invisible 
Sphere performs a loop in it and leaves in the direction of its entry (A), making the Sphere 
invisible. The singularity at the centre of the Invisible Sphere can be removed by expanding 
space around the centre within a given radius (B). 



This device has a spherically symmetric refractive index profile given by |4Tj 

n(r) = (q- 1/(3Q)^ if r < 1, and n(r) = 1 if r > 1, (3.1) 

where 

In the expression above we assume that the device has unit radius. It is easy to see that n(r) 
tends to infinity as r — > 0. This means that the speed of light goes to zero as it gets close to 
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the geometric centre of the device. This is extremely difficult to achieve in practice. 
Therefore, we make use of partial transmutation 
in 



39 



38 — a modification of the method presented 



to expand space with a non-constant factor in the immediate vicinity of the singularity 
within a given radius of the device (dashed circle in Fig. 4B). By this expansion we give more 
space for the propagation of light and thus ensure that the refractive index is finite everywhere 

see Appendix 



(i.e. light propagates with a non-zero velocity) within the entire device 38,39 



A.2) 



Maxwell's Fish Eye Lens 43,44 is an optical device in which all light trajectories trace out 



perfect circles (Fig. 5A). It has a spherically symmetric refractive index profile given by 



n(r) 



2n t 



(3.3) 



l + (r//) 2 ' 

In the expression above / is the radius of the "equator" of the device at which a circular mirror 
can be introduced (blue circle in Fig. 5B) such that all trajectories remain closed. 




Figure 5: Maxwell's Fish Eye Lens without and with the mirror. In Maxwell's Fish Eye 
profile all light trajectories are circular (A). When the Fish Eye Lens is surrounded by a 
circular mirror of appropriate radius, the trajectories still remain closed (B). 



Here ni is the value of the refractive index at the equator. It can assume any values that are 
greater than or equal to unity. A choice of large ni slows the light rays down, but it does not 
affect the trajectories overall, since the form of the trajectories in the device are dependent 
only on the relative distribution of the refractive indices. It is precisely this feature of the lens 
that we can exploit to avoid superluminal propagation in the inner branch of the cloak (see 
Appendix A.3). 
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4 Materials 



Now we turn to calculating the material properties of the cloak in physical space. We make 
use of the method described in |6,8,9|. We utilise the ideas of transformation optics [2-10 



which establish a connection between spatial geometries and dielectric media. We calculate 



the permittivity and permeability tensors in bipolar coordinates (see Appendix A.5). The 
tensors obtained have the following form: 

t\ = A = diag(e CT , e T , e z ) . (4.4) 

These two tensors completely describe the electromagnetic response of a dielectric material 
and their eigenvalues (i.e. e a , e T , and e z ) can be engineered into metamaterials to implement 
the cloak. The tensor eigenvalues and the refractive index values are related as [8l|9] 



n a = y/e T e z , n T = y/e a e z , n z = y/e a e T . (4.5) 

It is clear from the expressions above that the refractive index values could be greater than 
one even if not all of the individual tensor eigenvalues are greater than unity. However, in 
order to avoid anomalous dispersion and absorption at all frequencies, it is important that 
the individual tensor eigenvalues are also greater than unity. 

Through the appropriate positioning of the building blocks of our electromagnetic geometry 
and the adjustment of their parameters, all tensor eigenvalues can be brought above unity. 
We calculated the distribution of the tensor eigenvalues e a ,e T , and e z for the geometry shown 



in Fig. 3. The detailed calculations can be found in Appendix A.5 We plotted the eigenvalue 
distributions in Fig. 6. 




Figure 6: Distribution of the tensor eigenvalues within the combined cloak. The boundaries of 
the stretched region of the cloak are shown in black within the partially transmuted Invisible 
Sphere. All tensor eigenvalues are greater than or equal to unity and finite. 
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We found that e a G [1,9.428], e T G [1,73.21], and e z G [1,2012.12]. It is clear that all values 
of the tensor components are now greater than or equal to one, and finite. Therefore, all 
refractive index components are also greater than or equal to unity [see equation (4.5) . 
We note that some of these values are quite extreme and might be very difficult to achieve 
in practice. However, this cloak is just an example, which demonstrates that full invisibility 
cloaking is possible without superluminal velocities. The range of material parameters could 
be optimized by changing the topology of the cloak slightly (see Appendix A. 6) or by in- 
troducing further coordinate transformations to eliminate regions with high permittivity and 



permeability values (see 39 for an example). 

Figure 6 shows that the profile of the cloak is highly anisotropic within the inner branch of 
the cloak (demarcated by black lines) and isotropic outside it. In Fig. 7 we depict a light 
trajectory against the values of the e a component in the presence of the mirror. The light ray 
performs a loop within the transmuted centre of the Invisible Sphere, bounces off the mirror 
twice and leaves the device with its original direction restored. Objects placed in the white 
region behind the mirror are invisible. 

Finally, we note that even though this cloak is two-dimensional, the ideas behind it can be 
readily extended to construct three-dimensional version of it. A three-dimensional version 



of the non-singular expansion of space is given in 32 , and Maxwell's Fish Eye as well as 



the partially transmuted Invisible Sphere are both three-dimensional profiles. Naturally, the 
parameter ranges will be different for the three-dimensional cloak, but through appropriate 
optimization superluminal propagation can be avoided in the entire cloak. 



5 Conclusion 

In conclusion, we have combined the partially transmuted Invisible Sphere with Maxwell's 
Fish Eye Lens to obtain a cloak whose permittivity and permeability tensor eigenvalues are 
everywhere greater than or equal to unity and finite. This example demonstrates that full 
invisibility cloaking is possible without superluminal propagation and anomalous material 
requirements. 
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Figure 7: A light trajectory is shown against the distribution of the e a values. The light ray 
enters the device, completes a loop, bounces off the mirror twice and leaves the cloak with its 
original direction restored (A). Figure B gives a closer view of the vicinity of the inner branch 
of the cloak. Objects placed within the white region are invisible. 
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A Appendix 



A.l Branches in bipolar coordinates 

We use bipolar coordinates and the expansion of space to establish the two virtual sheets 
depicted in Fig. 2. We parameterise physical space by bipolar coordinates, which can be 
related to the Cartesian coordinates by 

a sinh r a sin a 

x = — , y= — . (A.6) 

cosh r — cos a cosh r — cos a 

In the expression above a is the radius of the boundary of the stretched region of the cloak 

(i.e. half the length of the branchcut). Working in the complex plane we can rewrite these 

relations succinctly (equations (S29 and S30) in |33|) 

fa + ir\ 

x + ly = 1 a cot I — - — J, (A. 7) 

with the inverse relation 

a + i r = 2 i arcoth ( x + iy \ (A.8) 



a a 



V a 

We map the physical plane (Fig. 2A) onto two virtual sheets (Fig. 2B) through the coordinate 
transformation given by (equation (S34) in [33] ) 

a for \a\ < — 

[Aa 2 , , \ c 7T , , (A.9) 

sgn a 3 cr + tt ) for — < \a\ < n . 

Here a' and r are the new set of coordinates parameterising virtual space. While a runs from 
—7i to tt, covering the plane once, <j'{<j) runs from — 2tt to 2tt and covers the entire plane twice 
thus giving rise to two virtual planes (Fig. 2B). Outside the black circle in Fig. 2A and 2B 
(corresponding to |a| < 7r/2) a and a' agree since that region undergoes no transformation. 
The two branches of virtual space are connected by the branchcut (yellow line in Fig. 2B), 
through which light can pass from one sheet onto another. In virtual space the bipolar and 
Cartesian coordinates are related by the expression 



a' + ir' = 2iarcoth( V ^ 1 ^ . (A. 10) 

It is useful to note that the expansion of space can be reversed using the inverse relation 
between a and a' (S35 in [33]): 

a' for \a'\<- 

vW) = {sma'{„ . ,. , A . tt , ! _ (A.ll) 



^3vr + ^1671- 1 <r' | - 7tt 2 ^ for ^ < \a'\ < 2n 
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A. 2 Invisible Sphere 



The refractive index profile of the Invisible Sphere (Fig. 4A) is given by equations (3.1) and 



( 3.2[ ) |4l] . In these expressions we assume that the device has unit radius. It can be shown 



that n(r) tends to infinity as r — > and diverges as ~ r ~ 2 / 3 near the geometric centre of the 
sphere. 

We remove the singularity at the centre by the method of partial transmutation |38|. We 
expand space around the singularity within a finite radius b to give enough space for the 
propagation of light. Thus we bring all refractive index values within a finite range. Our 
choice of b for the transmutation will depend on the particular geometrical set-up used (see 
section A. 6 on optimization), but as 38 shows, b will always have to be smaller than 0.2489 
to avoid tensor eigenvalues that are smaller than one. 

Since n{r) of the Invisible Sphere diverges as ~ r - 2 / 3 ) partial transmutation requires that we 
apply the following transformation to the radial coordinate: 

R{ r ) = — — r p+1 if r ^ b, and R{r) =r if r > b . (A. 12) 

In the equation above M is the normalisation constant that ensures that the transformation is 
continuous at the boundary of the transmutation r = b. Therefore, it follows that J\f = 6~ 2 / 3 . 
The material properties of the partially transmuted Invisible Sphere (Fig. 4B) - that would 
mimic the coordinate transformation described above - can be calculated in terms of spherical 
coordinates from the formula derived in [39] (see also |38|): 

• ■ . /nr 2 dR ndr ndr\ . . 

e\ = \i\ = &a%{e R ,e B ,e£ =diag(^— — ,— ,— J. (A.13) 

In the above expression n is the refractive index profile of the Invisible Sphere and R is given 



by equation (A. 12). 



For the particular partially transmuted Invisible Sphere positioned on the upper sheet of the 
virtual geometry in Fig. 3B, we chose the radius of the transmutation to be b = 0.075 to 



keep all combined tensor eigenvalues of the cloak greater than one (see section A. 6 on opti- 
mization). We plotted the transmuted values of £r and Eq = e^ in Fig. 8. Clearly, all values 
are finite and greater than one. Note that the region within the radius of the transmutation 
is anisotropic, whereas the outer region has been unaffected by the transmutation and thus 
remains isotropic. 

Also note that there is an alternative way of transmuting singularities without inducing anoma- 



lous material properties 40 . However, this method reduces the refractive index in the entire 



profile of the transmuted lens and, therefore, it is unsuitable for our purposes. 
A. 3 Maxwell's Fish Eye and the mirror 

The refractive index profile of Maxwell's Fish Eye Lens (43l|44] (Fig. 5A) is given by equation 



(3.3). We can solve Hamilton's equations numerically (see section |A^4 on Hamilton's equa- 
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Figure 8: Values of the radial and angular tensor components in the partially transmuted 
Invisible Sphere, where the radius of transmutation is b = 0.075. All eigenvalues are greater 
than one and finite. The index profile of the Invisible Sphere outside the radius of the trans- 
mutation remains unaffected. 



tions) to obtain the trajectories depicted in Fig. 5A. 

Maxwell's Fish Eye can be surrounded with a circular mirror of radius I (where I is the radius 
of the 'equator' of the device) such that all trajectories remain closed (Fig. 5B). To model the 
reflection of rays by the mirror, we express the light trajectories in the x-y plane in terms of 
the complex coordinate z'(t) = x'(t) + iy'{t) and apply the transformation 33 



I 



for 



> /. 



(A.14) 



{z'/iy 

In the expression above * denotes complex conjugation. 

For the particular Fish Eye Lens positioned on the lower sheet of the virtual geometry (Fig. 3B) 
/ = 2a, where a is half the length of the branchcut (see sections A.l and A. 6), and ni was 
chosen to equal 5 to make sure that all tensor eigenvalues are greater than unity within the 
inner region of the cloak. The Fish Eye is centered around (x',y') = (— a, 0). This point 
corresponds to one end of the branchcut. 



A. 4 Hamilton's equations 

In this section we describe how to obtain the light trajectories in physical space. To get 
the trajectories in the outer branch of the cloak (blue and black region in Fig. 3A), we first 
solve Hamilton's equations in the Invisible Sphere (Fig. 4A). We mimic the effect of the 
partial transmutation by transforming the radial components of the trajectories according to 



equation (A. 12) (Fig. 4B). Finally, we open up space within the Sphere along the branchcut 



using equation (A. 9), which cuts up the trajectories into two separate parts and pushes them 
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apart while the red region gets inserted in between them (Fig. 3A). The two separate parts of 
the trajectory in the blue region will be joined by the trajectory in the inner cloak. 
To obtain the trajectories in the inner cloak (red region in Fig. 3A), we solve Hamilton's 
equations in Maxwell's Fisheye profile with I = 2a and rii = 5. The initial conditions for the 
solution of Hamilton's equations are provided by the trajectories in the outer cloak, where 
the incident wavevector has to be modified to take refraction into account at the branchcut 
(see section A.6). Next we shift the trajectories such that they are centred around one of the 
endpoints of the branchcut (x', y') = (—a, 0). We switch from Cartesian to bipolar coordinates 



using equation (A. 10). Then we shift the values of a' from the intervals [— 7T, 0] and [0, n] to 
[7r,27r] and [— 2n, —it] respectively and reverse the expansion of space using equation ( |A.ll ). 
Thus we obtain the connected physical trajectories depicted in Fig. 3A and Fig. 7. 
We find it convenient for ray tracing to rewrite Hamilton's equations in a complex form. Using 
the dispersion relation u> = ck/n and setting c = 1, Hamilton's equations can be expressed as 
(equation (4.10) in (9|) 

dr lk 

dt n \ k\ ' 



(A.15) 



and 



dk 

~dt 



\k\ 



n- 



Vn , where Vn 



dn 



d\r\ 



We rewrite equation (A. 16) as 



dfc 

dt 



\ k \s(\ r \)n with s(l r l) 



1 dn 
n 2 d\r\ 



(A.16) 



(A.17) 



We solve equations (A.15) and (A.17) in terms of the complex position and momentum coor- 
dinates 

z(t) = x(t) + iy(t) and k(t) = k x (t) + ik y {t) . (A.18) 



Therefore, we substitute z(t) for r and kit) for k in equantions ( A.15[ ) and (A.17) and thus 
we obtain the complex form of Hamilton's equations 



dz 
dt 



k 



dk 
dt 



n(\z\)\k\> 

These equations are numerically stable. 



1*1 g(M) 



with g(|z|) 



1 dn 

n 2 d\z\ 



(A.19) 



A. 5 Calculation of material properties 

The material properties of the cloak can be calculated using the recipe described in [6|[8j[9]. 
We proceed in the following steps: (i) we write down the line element ds that describes our 
virtual geometry, (ii) express ds in terms of the coordinates that parameterise physical space 
(a, r, and z), and (iii) read off the metric tensor of the geometry from ds. Finally, (iv) we 
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calculate the permittivity and permeability eigenvalues from g^ using equations (8) and (9) 
from [6): 

v 7 ^ 



(A.20) 



In the equation above g 1 ^ is the inverse of gij, g is the determinant of gij, is the metric of 
the background geometry (here: empty space described by bipolar coordinates) and 7 is the 
determinant of 7y. In all our calculations we assume that the device is impedance matched 



to vacuum and so e\ = n\ [61. 



The cloak has two different branches that undergo different transformations, therefore, we will 
treat them separately. 



A. 5.1 Outer branch 

First, we calculate the material properties of the outer branch of the cloak (black and blue 
region in Fig. 3A), where tc/2 < \a\ < 3n/A. This physical region corresponds to the upper 
sheet of the virtual geometry (Fig. 3B). The upper sheet of the device carries the refractive 
index profile of the Invisible Sphere n is . Therefore, the line element in spherical coordinates 
is given by: 

ds 2 = n 2 s (dr 2 + r 2 d# 2 + r 2 sin 2 # d0 2 ) . (A.21) 

We perform the radial expansion given by equation ( A.12[ ), which corresponds to partial 
transmutation 



ds 2 



dr 
dR 



dR 2 + r 2 d6 2 + r 2 sm 2 6 d(f) 2 



(A.22) 



Next we change from spherical to Cartesian coordinates 

2 



ds 2 



n r . 



+ 



+ 



dR) \dx) 
dR) \dz ) 



de 

dx 

de 

dy 

de 

d~z 



+ r 2 sin 2 ^ 



+ r 2 sin 2 # 



+ r 2 sin 2 ^ 



dx 



dy 

d(j) 
dz~ 



dx 2 
dy 2 
dz 2 



(A.23) 



We define the expressions in square brackets in front of dx 2 , dy 2 , and dz 2 as e x , e y , and e 2 
respectively and obtain 

ds 2 = n 2 s (e x dx 2 + e y dy 2 + e z dz 2 ) . (A.24) 
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We switch from Cartesian to bipolar coordinates using equation (A.6) 

2 



ds = nz 



dx 
da' 



dy_ 
da' 



fdz 



( dx\ ( dy\ ( dz 



da 



dr 2 + e z dz z 



(A.25) 



Again, we define the expressions in square brackets in front of da' 2 and dr 2 as e a > and e T 
respectively and obtain 

ds 2 = n 2 s (e a > da' 2 + e T dr 2 + e z dz 2 ) . (A.26) 
Finally, we change from a' to a 

2 



ds 2 = n 



e ^'(^p) da 2 + e T dr 2 + e z dz 2 



(A.27) 



da' 



We define e CT = e a i — and deduce the metric tensor from the above line element 
\daj 

9ij = n 2 s diag(e CT , e r , e z ) . 
The determinant of this tensor is g = nf s e a e T e z and its inverse is given by 



t'J 



1 A ■ 1 1 1 

— diag — , — , - 



IK 



£<T ^-Z 



The metric tensor of the background geometry is given by 

7 ifc = diag(7 CT ,7 T ,7 2 ) . 
The determinant of 7^ is 7 = ^alrlz, where 



la = 7 



(cosh r — cos a) 2 



and 7 Z = 1 . 



(A.28) 



(A.29) 



(A.30) 



(A.31) 



In the expression above a is the radius of the stretched region of the cloak [see equation (A.6) 



We can now substitute the relevant expressions into equation (A. 20) and obtain the permit- 
tivity and permeability tensors 



H\ = diag(e (T , e T , e z ) = n IS diag 



7(T^T^2 / ^ct7t^2 



^alrlz V Icy^-rlz V T<t7t^z 



totrlz 



(A.32) 



The permittivity and permeability tensors given above describe the electromagnetic properties 
of the outer branch of the device. Note that we expressed the material properties in terms 
of the mixed (co- and contravariant) tensors, since they are directly related to the refractive 
indices [see equation ( |4.5 )] and thus have physical significance 39 . 
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A. 5. 2 Inner branch 

Next, we calculate the material properties of the inner branch of the cloak (red region in 
Fig. 3A), where 37r/4 < \a\ < n. This physical region corresponds to the lower sheet of the 
virtual geometry (Fig. 3B). The lower sheet carries the refractive index profile of Maxwell's 
Fish Eye n M given by 

5 • 2 

nM(X '' ^ = l + {(x' + a)/(2a)y + {y>/(2a)r ^ 
Therefore, the line element in bipolar coordinates is given by 

ds2 = ""{( (co S hx-co S ^ )^ + dT2 » +dZ 1- < A ' 34 > 
To simplify the line element we define A s = A T = (a 2 /(coshr — cosa) 2 ) and obtain 

2 



ds 2 = <| \„da 2 + A T dr 2 + dz 2 j . (A.35) 

Next we bring the values of the sigma coordinate from the ranges [—it, 0] and [0, it] to [it, 2tt\ 
and [— 2n, —it] respectively using the function 

//-\ [o — 2it for a > /A „„ N 

a'(a) = < _ _ , _ ^ _ (A.36) 

v ; \ a + 2 vr for cx < . v ; 

We obtain 

ds 2 = ni^^^da' 2 + A T dr 2 + dz 2 } = n 2 M ^X,da' 2 + A r dr 2 + dz 2 } , (A.37) 
where the last equality follows from the fact that (da /da') 2 = 1. Finally, we switch back from 



a' to a to reverse the expansion of space using equation (A. 9) and obtain 



Off 



ds 2 = n 2 M \\^[ — ) da 2 + X T dr 2 + dz 2 \ = n 2 M \ \ a da 2 + X T dr 2 + X z dz 2 } . (A.38) 



In the expression above the last equality follows from defining the factor in front of da 2 as A CT 
and defining X z = 1. The metric tensor of the background geometry is given by 

K jk = diag(K CT) k t , K g ) . (A.39) 

The determinant of where 

a 2 

— K T — — and k z — 1 . (A. 40) 

(coshr — cos a) 2 
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We substitute the relevant expressions into equation (A.20) and thus obtain 



e ki ^ = f* k = diag(s ff) s T ,s,) = ^dxag^^^,^^-, ^ j . (A.41) 

Fig. 6 shows the distribution of the tensor eigenvalues within both the inner and outer branches 
of the combined cloak. For the parameter ranges we get e a G [1,9.428], e T G [1,73.21], and 
e z G [1,2022.12]. 



A. 6 Optimization 

In this section we determine the optimal parameters for the topology that combines the 
partially transmuted Invisible Sphere and Maxwell's Fish Eye Lens (Fig. 3B). We emphasize 
that this geometry is merely an example, which demonstrates that full invisibility cloaking 
is possible without superluminal velocities. This geometry could be modified in various ways 
(e.g. by changing the topology slightly or by introducing further coordinate transformations) 
to find a set-up with less demanding material parameters. 

While positioning the branchcut (yellow line in Fig. 9) within the Invisible Sphere (purple circle 
in Fig. 9), we have to ensure that all index values are raised above one within the stretched 
region of the cloak and that total internal reflection is avoided within the cloak. Without the 
background of the Invisible Sphere, the lowest eigenvalue of the stretched region of the cloak 
(blue region in Fig. 3A) is 0.333 in the e a profile. This value lines the entire length of the 
branchcut. Therefore, the branchcut has to be placed in the region where n IS > 3 (within the 
green circle in Fig. 9). 

Total internal reflection may occur when a light ray passes from an optically dense medium 
to a lighter one. To avoid total internal reflection, we have to ensure that the refractive index 
along the branchcut is higher in the inner cloak than in the outer cloak. Since the expansion 
of space was achieved by a smooth function, the expansion by itself would give rise to a 
smooth index profile. Therefore, total internal reflection will depend solely on the values that 
the Invisible Sphere and Maxwell's Fish Eye assume on the two sides of the branchcut. The 
Fish Eye profile changes between 5 and 10 and reaches a value of about 8 at the midpoint 
of the branchcut. Calculations show that total internal reflection is avoided and the size of 
the branchcut is maximised when the branchcut is oriented as shown and is placed within the 
blue annulus in Fig. 9, where the red circle seals off the region where n I3 > 7.89. Assuming 
that the Invisible Sphere has unit radius, we obtain n = 0.08 and r2 = 0.288 for the radii of 
the red and green circles respectively. This yields 0.554 for the length of the branchcut (i.e. 
0.277 in equations flA~6| \K1\ |A~tj |A~3~T] and [AT40| . 



a 



Part of the stretched region (black circle in Fig. 9) will now lie outside the blue annulus. 
However, calculations show that e a goes sufficiently quickly to unity towards the boundary of 
the stretched region to ensure that the combined index values do not dip below one. 
Finally, we have to choose a sufficiently small radius for the transmutation of the Invisible 
Sphere to raise the combined eigenvalues within the transmuted region above one. We found 
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Figure 9: The position and size of the branchcut in the Invisible Sphere. The blue annulus 
highlights the region where the branchcut can be placed such that all tensor eigenvalues 
are greater than unity and total internal reflection is avoided. The yellow line indicates the 
orientation and maximal size of the branchcut. The black dashed circle shows the boundary 
of the transmutation of the Invisible Sphere. 



that b = 0.075 is the greatest and thus optimal transmutation radius that achieves this. 
Note that refraction will occur at the branchcut when light crosses from the outer to the inner 
branch due to the differing values of the Invisible Sphere and the Fish Eye profile on the two 
sides of the branchcut. However, due to the particular symmetry of the device, when the light 
ray crosses the branchcut for the second time, the effect of the refraction is exactly reversed 
and thus the trajectory is not distorted overall. 

We note that, for the purposes of optimization, the refractive index of the inner cloak of our 
device can be replaced by an index profile that is identical to that of the inner cloak of the 
Non-Euclidean cloak described in 32 and 37 except for an overall multiplicative constant of 
5.274. This constant has to be introduced to ensure that all index values in the inner cloak 
are above unity. Following optimization, we find that the eigenvalues for this construction are 
E(j G [1, 11.65], e T G [1, 59.2], and e z G [1, 944.9], which shows that the range of the parameters 
is already a factor two smaller than the range for the cloak with the Fish Eye profile. However, 
the size of the invisible region also becomes smaller for this cloak. It is generally true that 
the greater the invisible region, the higher tensor eigenvalues are needed. 
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